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An important open problem in quantum information theory is the question of the existence of 
NPT bound entanglement. In the past years, little progress has been made, mainly because of the 
lack of mathematical tools to address the problem, (i) In an attempt to overcome this, we show 
how the distillability problem can be reformulated as a special instance of the separability problem, 
for which a large number of tools and techniques are available, (ii) Building up to this we also 
show how the problem can be formulated as a Schmidt number problem, (iii) A numerical method 
for detecting distillability is presented and strong evidence is given that all 1-copy undistillable 
Werner states are also 4-copy undistillable. (iv) The same method is used to estimate the volume 
^— ^ ' of distillable states, and the results suggest that bound entanglement is primarily a phenomenon 

^SJ ' found in low dimensional quantum systems, (v) Finally, a set of one parameter states is presented 

1— ( , which we conjecture to exhibit all forms of distillability. 

a^ : 
^ . 

I. INTRODUCTION 

cn ■ 

>'■ 

lO ■ Only recently the mathematical definition of entanglement in quantum information was rigorously supplemented 
by a physical interpretation. The definition, as introduced by Werner is well known. A bipartite system p G 



o 



C{Ti.A ^Ti-B) is called separable if and only if p can be expanded as 

O • with Pi > 0. If this is not possible, p is called entangled. When a state cannot be written in a separable form, does 
I ' . . that mean it cannot be constructed locally? For a single copy of a state this is obvious P|, while in the asymptotic 
^ , regime this was only recently proved 0]. 

^ ' Considerable effort has been devoted in quantifying the amount of entanglement present in a state, usually with 
Cr' the aid of entanglement measures M The principal physical demands on such an entanglement measure is that 
^ it vanishes on separable states, and that it is non-increasing under local operations and classical communication 
I (LOCC). However, it is by no means evident what 'entanglement present in a state' means. There are two obvious 
^ , physical options to this: we can either mean the amount of entanglement used to construct the state or the amount 
■ ■ ■ ' of entanglement we can recover from the state. Usually this is done in the asymptotic regime and entanglement is 
measured with reference to the singlet state \^) = :^(|00)-|-|11)). The two associated entanglement measures are then 
called the entanglement cost Eq and the distillable entanglement Er). A celebrated result of quantum information is 
that both values coincide for pure states and are equal to the von Neumann entropy of the reduced density operator 
Q. For mixed states it was shown that any entanglement measure E should satisfy E]j < E < Ec 0. The results 
of Ref. 01 we mentioned earlier implies that Ec > for all entangled states. The question whether Ed > for all 
entangled states was answered negatively. Indeed, Horodecki et. al. showed that there exist entangled states from 
which no entanglement can be distilled at all When Ed > we call the state distillable, otherwise it is called a 
bound entangled state. 

The question was then to classify all bound entangled states. This is also known as 'the distillability problem' and 
is the main focus of this paper (see problem 2 in Q). We start off by reviewing what is known. The next theorem is 
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crucial. 

Theorem I.l (Horodecki et al. (i) All entangled two qubit states are distillable. 

(a) An arbitrary bipartite state p acting onT-Cy^^T-Cg is distillable if and only if there exist projectors P : 7i®" 7^2 
and Q ; 7i^" — ^ TI2 and a number n, such that the state 

p' ^{P(g>Q)p®"{P(g>Q)^ (1) 

is entangled. Since p' is a state acting on 7^2 ® 'H2, this means that p' needs to have a negative partial transposition: 
{p'f^ < 0. 

If such an n exists we call p pseudo-n-copy distillable or in short 71-distillable. The prefix pseudo reflects the fact 
that if we project upon such a subspace we are only half way through our distillation process. Indeed, in the next 
step we would like to repeat this procedure m times on batches of n copies of p, giving us m copies of the qubit pair 
p'. Finally we can use existing protocols to extract maximally entangled singlets from p'®™ |lOlllll| . 

The theorem is equivalent to saying that if a state p is distillable then we can find a Schmidt rank two vector 
i/j and a number n such that ("iAKp^lHJ^ IV") < 0- From this it follows that states with a positive partial transposition 
(PPT) can never be distilled (see [l3j for a direct proof). Usually the term 'bound entangled states' is therefore 
associated with entangled PPT states. The question remains whether there exist bound entangled states with a 
negative partial transposition (NPT). This problem can be reduced to the question whether all entangled Werner 



[isposition 

states are distillable |lLll3lll4l|. Recall that a Werner state acting on TCa ^Ti-B = (8) can be written as 

Here F — \ij){ii\ denotes the swap or fiip operator. These states are entangled when Tr[pwF) < or /? < 
— 1/d. A prominent property of these states is that they are the only states which are invariant under local unitary 
transformations of the form U ®U . Hence, any state p can be transformed into a Werner state by applying the so 
called twirl operation: 

pw^ I dUU®Up{U®U)\ 



where the integral is with respect to the Haar measure on U{d). Note that this transformation leaves the expectation 
value Ti{pwF) = Tr(pF) invariant. Now it is easy to prove that one can transform any NPT state p to an NPT state 
p' such that Tr(p'F) < 0. Applying the twirl then gives an entangled Werner state. 

The distillability of the Werner states has been studied in two papers [13. 0| . The authors were able to show that 
they are distillable when P < —1/2 and rt-copy undistillable for /3 > — + e„. Unfortunately, the range e„ goes 
to zero as n goes to infinity. It is however conjectured that the Werner states are undistillable for all f3 > —1/2. 
An important result in this context was obtained by Watrous jl^ who constructed a one parameter set of distillable 
states which are n-copy undistillable in some range. Supporting evidence in favour of the conjecture was provided 
[1^ in the form of numerical evidence for 2 and 3 copies for d ~ 3. Apart from intrinsic importance of the conjecture, 
affirmation would imply non-additivity and non-convexity of bipartite distillable entanglement T?] . The problem also 
has non-trivial consequences on the theory of positive maps |15|. 

A remarkable effect in the context of distillation, is activation of PPT bound entanglement. It has been proven 
that every state becomes 1-distillable by adding a PPT bound entangled state [liEi (see also 0). Conversely it 
has been shown that for every PPT state there exists a 1-undistillable state, such that taken together one obtains a 
1-distillable state 0. 

In the next section we will discuss the distillability properties of certain class of highly symmetric states, which 
includes the Watrous states and two copies of the Werner states. We will come back to this set of states repeatedly 
in the rest of the paper. Section HI is the main part of this paper, here we reformulate the distillability problem first 
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as a special instance of the Schmidt number problem, and using similar techniques as an instance of the separability 
problem. We suggest and discuss several approaches to tackle this specific separability problem. In section IV we 
outline our numerical method for detecting distillable states, we give a numerical estimate of the volume of distillable 
states for low dimensions and we provide strong evidence that all 1-undistillable Werner states are also 4-undistillable. 
In the Appendix we present a one parameter set of states which appears to exhibit all forms of distillability. In most 
of this paper we will omit the normalisation of density operators as they are not relevant to us. 



II. THE [/C/V'VF-INVARIANT STATES 

To start we will briefly recapitulate some properties of the so-called local symmetry groups and states invariant 
under such groups. For an excellent overview with plenty of examples, the reader is referred to Ref. [2^ . 

Let Go be a subgroup of the unitaries, or possibly the whole group of the unitaries. Then we can define the group G 
of the unitaries consisting of all pairs of the form U ®U' acting on a Hilbert space Ti = Ha ^Hb, where U G Gq and 
U' is some given function of U. The set of bipartite states left invariant by G is just the intersection of the state space 
with the commutant of the group G. Generally speaking, choosing Go sufRciently large, there will exist a finite basis 
of operators spanning the commutant. The set of bipartite states left invariant by the group G will also be denoted 
as f/C/'-invariant instead of G-invariant. An arbitrary state can projected onto an f/t/'-invariant state by twirling it 



'?g(p) = j dUU®U'p{U®U')\ 



here the integral is performed according to the Haar measure on Gq. In the introduction we have seen an example 
of such a local symmetric set of states, namely the Werner states, which are [/^/-invariant and spanned by 1 and F. 
The so called isotropic states 0,0 are [/J7*-invariant and a basis is given by 1 and P = l/^X^ij I**)OjI- 

From these two basic symmetry groups one can generate others by considering tensor products. Consider the case 
where we have two symmetry groups G — {U ® [/'} and K = {V V'} acting respectively on Hi = Ha ^ Hb 
and H2 = H'a ^ H'b- Let Bq and Bk be a basis for the UU'- and the V^F'-invariant states respectively. Then 
a basis for the t/J/'Fy'-invariant states acting on Hi <S) H2 will be given by Bq ^ Bk- In what follows we will 
number the systems belonging to party A with odd numbers and party B with even numbers. As an example, a 
basis for the UUVV-mva.riant states is given by the operators {I12 ® I34, ^12 ^ I34, 112 ^Fsa, F12 (8) -F34}. Imposing 
the extra condition that the coefficients of F-\2^ I34 and I12 (Xi-F'34 should be the same, we end up with the so-called 
UUVV F-mva,riant states as introduced in |22l |. where it was used as a counterexample to the additivity conjecture 
for the relative entropy of entanglement. 

The set of [/C/yVF-invariant states contains the Werner states and the Watrous states and hence is also ideal for 
the study of distillability. A convenient parametrisation is given by 

, l-2ed + 6d'^ 

/9 = I12 I34 + — 3— (I12 «) -F34 + ^12 I34) H -p^ -^12 ® -^34- 

d d^ 
The set of density operators is restricted by the following inequalities 

(d - if + 2ed(d + 6d^ > 0, 
_ 1 + 2ed - 5d^ > 0, 
{d + if - 2ed{d + 1) + 6d^ > 0. 

This set of states includes the Werner states (2 pairs in d (g) d for 5 = and 1 pair in d^ (g) d^ for e — 1/d) and 
the Watrous states (1 — 2ed + Sd^ = d^). The PPT states are just the separable states (see [i^). The states are 
entangled for e < or 5 < which follows from 



P^"" = Q12 Q34 + de{Pi2 (334 + Q12 «) P34) + 5d^Pi2 «) P34, (2) 
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with (3 = 1 -P. 

Let us now investigate the distillation properties of these states. To prove distillability, all we need to do is to find 
Schmidt rank two vectors ip such that (iplp^^lip) < 0. We present three such vectors, which we conjecture to detect 
all 1-distillable states: 

|^).l = |00)i2®(|00> + |ll))34, 

|^)b = |01)i2(»(|00) + |11))34, (3) 

Wc = J2 \^^)A\jj)B + + l))A\jij + 1))b. 

which detect distillability for the states respectively satisfying 

(f + U{ed - 1) + 2(1 - 2(d + 5d^) < 0, 
1 1 

'^d-r 

X 1 1 
'^d^-2- 

This set of states is shaded dark grey in Fig. 

Now let us derive results in the other direction, namely which states are undistillable? First let us show that we 
can find NPT states which are n-undistillable for arbitrary n. We will use the following inequalities derived in 



1 ^N-k 



(V'lg^^l^) > (1 



2\JV 



for any \tp) with Schmidt rank two. Now from jSJl it follows that all potential negative terms contain at least one P 
term. But we will also have a term (vE'jQ^^^I'J') > (1 ^ i)^^' this term can always dominate when we choose e 
and 6 small enough. Thus for each n, as long as we choose e or S small enough (but one or both negative), we obtain 
an n-undistillable NPT state. 

Next, let us derive states that are definitely 1-undistillable. From ^ and the given inequalities it follows straight- 
forwardly that the states satisfying 

(1 - 2/d)2 + min(4e, 0) + min(2(5, 0) > 0, 

are 1-undistillable. Note that this set does not touch the set of 1-distillable states. The reason for this is that although 
the inequalities are sharp, the sum of such inequalities are not. In some regions we can get a better bound by rewriting 
@ as 

p'^^ = li2 (g) I34 +{ed - l)(Pi2 (E) Q34 + Q12 «) P34) + (Sd^ - l)Pi2 «) P34, 

and we find that all states such that 

d^ + mm{4:d{ed - 1), 0) + mm{2{Sd^ - 1), 0) > 

are 1-undistillable. This second set of states is depicted in Fig. 

To conclude our discussion of the distillation properties we will discuss the Watrous states Q| of which the general 
form is given by 

p = li^2 ® 13,4 +^^^^(ll,2 (E)F3A + Fi^2 ® 13,4) + Fi^2 ® ^3,4, 
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FIG. 1: UUV V F-invariant states. All states satisfying e > and S > are separable. Shaded areas mark out distillable states 
(1,2 and oo- copy distillable). 



with l + l/d>e>l/d— 1. The states are entangled if and only if e < and definitely 1-distillable ii e < 1/d — 1/2. 
We will show that all entangled Watrous states are distillable. Suppose now we have two pairs of entangled Watrous 
states, the second pair having indices 5, 6 and 7, 8. Using the identities 

Tr((Fi,5®P2,6)(ll,2<»F5,6)) = l/d, 
Tr((Pi^5 P2,6)(Pl,2 ® l5,6)) = l/d, 
Tr((Pi,5 ® P2fi){Fi,2 ® ^5,6)) = 1, 

it can be verified that projecting upon Pi^s ® P2.6 will yield a new Watrous state p' with parameter 

/2(erf + d^_l)\ 

Since e' < e for 0>e>l/(i — 1, the state will be more entangled than the state we started from for each e < 0. One 
can repeat this protocol on many pairs until finally e < l/d — 1/2, at which point we obtain a 1-distillable state. More 
generally the protocol can be applied to the whole set of states; a straightforward but tedious calculation leads to 



e{dH + d^-l) 

5' ^ 



d^e'^ + d'^-l ' 
e2(d2_i) + rf2j2 



rf2g2 + ^2 - 1 ■ 

The states which are 2-distillable with this protocol are depicted in Fig. ^ Repeating the protocol recursively, it is 
not hard to show that all states satisfying 

3d2 + 4d-8 1 
^ 2d{d-2) ^+ d2' 

are distillable (see again Fig. In section IV. A we will give evidence that all the other states are probably NPT 
undistillable, but this of course this awaits an analytical proof. 



6 



III. A POSITIVE APPROACH TO THE DISTILLABILITY PROBLEM 

As mentioned in the introduction, the distiUabihty conjecture is equivalent to the statement that there exists no 
Schmidt rank two vector |?/;) such that 

for aU n > 1 and = 1— In other words, it seems hke affirmation of the conjecture would have to be in 
the form of an impossibility proof as opposed to a constructive proof. We will reformulate this conjecture in a more 
tractable form, namely as a special instance of the separability problem, for which a large number of tools are present. 
As a steppingstone we first show how to translate the distiUabihty problem into the problem of detecting Schmidt 
number 3. In the next subsection we then reformulate it as a separability problem. 



A. As a Schmidt number problem 

The Schmidt number of a quantum state has been introduced in as a generalisation of the Schmidt rank of 
a pure state. The Schmidt number of a state p is defined as the smallest number n such that p can be written as 
a convex combination of pure states with Schmidt rank n. Thus, separable states have Schmidt number one and 
entangled states have a Schmidt number larger than one. The problem of detecting the Schmidt number of a state 
has received little attention until now. A notable exception is Ref. |25| which developed the notion of Schmidt-number 
witnesses. When a state is 1-undistillable we have that Tr(|-0)(V'|p'^^) > for all Schmidt rank two states ip. Thus 
for every NPT 1-undistillable state p, is positive on Schmidt rank two states and is thus a Schmidt number 3 
witness. Analogously, 1-distillable states will give rise to Schmidt number 2 witnesses [2^. The symmetry of the 
Werner states allows for a dual approach, reformulating the problem of proving that a certain operator is a Schmidt 
number 2 (three) witness into the problem of detecting the Schmidt number 3 of a certain class of states. We will do 
so for one, two and n copies of the Werner states. 

Let us start with one copy of pw , for which the answer is known. As pw belongs to the set of C/?7- invariant states, 
p-^ will belong to the set of J7J7*-invariant states and hence will be invariant under the [/C/*-twirl: Tuu* {p^ ) — Pyy ■ 
From this follows that 

ii^lp^l^) = Tr{P^p^) = Tr{Tuu-.{Pi,)p^), 

here Tuu'{Pii) is the operator = after apphcation of the C/J7*-twirl. Thus we do not need to check over 

the whole set of Schmidt rank two vectors, but instead over the restricted set of [/[/*-invariant states with Schmidt 
number 2. The Schmidt number of the C/J7*-invariant states or the isotropic states is well known 13 E^- If we 

parametrise the isotropic states as = 1 +aP then pa has Schmidt number k when 

dikd-l) 

From this follows that pw — 1 is one-distillable if and only if Tr(p^pQ) < 0, with a — ^^^§35^- Going through 
the algebra we recover that all Werner states pp with (3 < —1/2 are 1-distillable. 

Let us now look at two copies of the Werner states p^^; as pointed out before, these states belong to the larger 
class of the [/[/y^i^-invariant states. Thus the relevant dual set is the set of the UU*VV* F-\nvaiia.\A states. For 
convenience we write the operators in the order of the indices 1,2,3,4 and omit these indices. With this in mind we 
can parametrise the UU*VV* F-\nvaiiaiA states as 



p = Q®Q + x{Q ® P + P ®Q) + yP ® P, 
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with x,y > 0. These states are separable for cP — 2d{x — 1) + (1 — 2a; + y)<0 and cP — (1 — 2x + y) < as depicted 
in Fig. 12 The set of Schmidt number 2 states contains at least the convex hull of the points 

/ (3d-4)(d+l) 2{d+l)\d-l) \ 
\ 2d - 4 ' d-2 J 

-(^•°) 
-(°.^) 

obtained by twirling the Schmidt rank two vectors from equation Q. All the states lying above the C^-line have 
a Schmidt number larger than two. This follows from the fact (see section II) that the operator ^P(8)P is 

positive on Schmidt number 2 states. If the Werner states are 2-undistillable then also all the states lying to the right 
of the ^i?-line have Schmidt number at least three. This follows easily by evaluating the expectation value of p at 
(1 — It is important to note that it is sufficient to show that all the states p on for instance the i?_D-segment 

(see Fig. ^ have Schmidt number 3. This can be seen as follows: let p approach E, and suppose we can show that 
each p has Schmidt number 3. This would imply the existence of a hyperplane W separating p from the set of Schmidt 
number 2 states. For p arbitrarily close to E this hyperplane would be parallel to the AB-segment, otherwise cutting 
it. Therefore this would show that all the states lying to the right of the Ai3-hyperplane have Schmidt number at 
least three. 

For the general case we need to consider the set of UiUi ■ ■ ■ J7„[/*F-invariant states. Here the subindex in Ui refers 
to the subsystem the unitary operator is acting and the F denotes any permutation of the subsystems. In what follows 
we will call these states C/it/*F- invariant states. We parametrise them as 

p = Q®" + xi {p (g> Q®"-i + Q(g>P(g> g®"-2 + •••)+ 

cc2(P®2 (8) Q®"-2 + P®Q®P® Q®"-' + •••) + •••+ XnP®", 
where we have found it convenient to now use the normalised Q — Q/(d^ ~ !)• The relevant hyperplane in this case 
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is given by 

^•[(Q.(ij)p)%]^i.i:(';)(-|)'^.^o. 

We can also write this as 1 + X]"=i = ^ with = (| — l)'^;^. From now on we wiU continue to work in 

these normahsed variables. Next we will generalise the idea developed for two copies. First we need to show that the 
hyperplane is spanned by Schmidt number 2 states. Then in order to check distillability, it will be enough to find 
the boundary between two and three Schmidt number along a line from the origin to an interior point of the points 
spanning the hyperplane. 

An independent set of Schmidt number 2 states spanning the hyperplane is easily obtained as follows (compare to 
the case for two copies). Let \i>)\ — -^lOl)*^""-^® (|00) + |11)). Twirling this state will yield an C/^ [/* F-invariant state 
with coordinates (i,0, •• • ,0). Similarly, twirling 1-0)2 = ^|00)(K)|01)®""2^(|00) + |11)) wiU yield an [/,[/* i^-invariant 
state with coordinates (x\ 7^ 0, i2 7^ 0, 0, • • • ,0). In general 

|V^)fe = ^100)^^= ® |oi)®"-fc-i ^ (|oo) + 111)) 
v2 

will yield a state with coordinates xi ^ for i < fc + 1 = for « > /c + 1. It is evident that all points will lie on 
the hyperplane and that they form an independent set, spanning the hyperplane. An interior point can for instance 
be obtained from the first point, as (i + e(l — i), e, • • • , e) for sufficiently small e. It can readily be verified that this 
point belongs to the hyperplane by using the identity Y^l^^ = " ^ 1- 

Thus n-undistillability of the Werner states beyond the 1-distillability boundary is equivalent to the statement that 
the I7il7*i^-invariant state with coordinates + e(l — i) + 5, e, • • • , e) has Schmidt number 3 for e > small enough 
and aU (5 > 0. 



B. As a separability problem 

This section contains the main result of this paper, namely the casting of the distillability problem as a special 
instance of the separability problem. An important tool in this, is the following result. 



Theorem III.l (Kraus, Lewenstein and Cirac |26|). Let Pi he the projector onto a maximally entangled state 
acting on Hi = Ti.Ai = (8" C^. Then for an arbitrary operator X acting on Tii = 'Ha2 ®'Hb2 we can define 

Wx^P2<E)X^^. (4) 

A state p acting on Ti.2 is n-undistillable if and only if Wpfsn is an entanglement witness. 

This theorem can readily be seen from the following lemma 

Lemma III.l. Let a be a positive operator with Schmidt number N > I acting on Tii — 'S) Ti-Bi '^'^'^ V '^'^ 
operator acting on Ti.2 = 'HA2 ® ^-82 positive on states with Schmidt number KN . Then the operator a ®rj acting on 
Til ® Ti.2 — Ti.A ® Tin is positive on states with Schmidt number K . 

Proof. It is clear that it is sufficient to prove the lemma for pure states a = |0)(0|. So let 

N 



and take an arbitrary Schmidt rank K state 



= '^(i}i\ai)A^\h)B^, 



K 

^^j\(^j)A\fj)B. 
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FIG. 3: The state space of the ?7?7* VV^*-mvariant states. The point C hes in the xz-plane such that the points A, B, C and E 
lie in the same plane. 



Then we need to prove that Trd^p) {ip\a (g) 77) > 0. This trace operation can be composed of tracing out the first pair, 
and then the second, as Tr(-) = Tr2(Tri(-)). Then making use of the identity Tri(C(Ai (g) B2)) = Tri(C(Ai (g) l2))B2 
we have 

Tr(|V')(V'k ^V)^ Tr2(Tri(|V')(^|(a (g la))?]) 

Now Tri(|?/;)(V'|(crg)l2)) is the projector onto the pure Schmidt rank iiTiV state 17)2 — (^IV') —^ij'l>i'^j{o-i\ej){bi\fj) £ 
7^2 from which the lemma follows. □ 



Special cases of this lemma have appeared in the literature over the years jl2l Il8l 1261 127| . Next we will apply 
Theorem IIH . 1 1 to the Werner states for one and n copies and making use again of the local symmetry, we will present 
a dual positive formulation of the conjecture. The analysis will be very analogous to the reformulation as a Schmidt 
number problem and a continuous comparison of this section with the previous is very useful. 

For one pair, we need to prove that P2 (g (1— |P) is an entanglement witness. As before, it will be sufficient 
to characterize the subset of the separable states of the general UU*VV* invariant states. Here U acts on a two- 
dimensional Hilbert space. We parametrize the C/[/*yF*-invariant states as 

P^Q2<^Q + xP2(dQ + yQ2 ®P + zP2®P, 

with x,y,z > 0. The separable states are a subset of the states with positive partial transpose, which satisfy the 
inequalities 

ABCE-plane 
ECD-plane 
EAD-pla,ne. 





- ix 
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— X 






+ 1 


+ y 
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- 1 






- 1 


+ y 



Thus the states with positive partial transition are contained in the polyhedron spanned by the points OABCDE as 
in FigOl By twirling the pure separable states |0101), |0100), ^|^-^g |0001) and |0000) one obtains the states 
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represented respectively by the points 0,A,C,D and E. Note that B is not in this list. Indeed, as we know that 
is an entanglement witness, all states satisfying z > 2x/(d — 2) are entangled. In particular the states in 
the tetrahedron spanned by the points ABCO are PPT entangled. Conversely, knowing that the polyhedron ABCO 
is PPT entangled immediately proves that P2 <S) (1 — 5^) is an entanglement witness. 

For two copies, it is sufficient to study the set of C/C/*(ViF]*V2V^*P)-invariant states parametrised as 

P = Q2® {Q ® Q + yi{P ® Q + Q ® P) + 2/2 -P ®P) + P2® {xoQ Q + xi{P Q + Q P) + X2P «> P) 

We will not attempt to completely classify the separable subset, instead it is enough to look at what happens in the 
neighbourhood of the hyperplane 



Tr 



P2 (Q+ (1 - ^)P) p]=Xo- {d-2)xi + (^) X2 = 



0. 



We now show that this hyperplane is spanned by separable states. Consider the following pure separable states and 
their coordinates {yi,y2,xo,Xi,X2) after action of the UU* {V1V1V2V2 F) twirl: 

|01)|01)|01): (0,0,0,0,0) 

|01)|00)|01) : |^^^^i-^,0,0,0,oj 

101)100)100): [-^^,,^^,M0^ 



i,j=0 

1 



2(rf-l)' ' 3d ' 6d{d-l)' 3d{d-l) 



^ |ii)|ii)|01) : 0,0,^,1,0 . 



i,i=o 



These coordinates can be verified after a tedious but straightforward calculation. 

A point in the interior of the convex hull of these points can be obtained by averaging these coordinates. In this 
case one obtains the state with coordinates xq = 2{d — 2)/(15rf), xi = {5d — 6)/{30d{d — 1)), X2 = 2/(15d(d — 1)), 

yi = d{2d - 3)/(5(d - l){d'^ - d - I)) and 7/2 = l/(5(d^ - d- 1)). From this follows that the Werner states are 2- 
undistillable if and only if the UU* {V1V1V2V2 F)-mva,na,n.t states with coordinates {yi,y2,Xo, Xi+e, X2) are entangled 
for all e > 0. 

Let us now move to n copies. The relevant set of states is the set of ?7J7*(ViV^*F)-invariant states 

P = Q2 ® Q®" + Q2 ® [yi{P ® Q®""' + Q®P® 0®"-' + •••)+ i/2(P®^ ® Q®""' + •••) + ••• + 2/„^'®""']+ 

P2 ® [a^oQ®" + xi{P ® 0®"-i +Q®P® Q®"-2 + ...)+...+ a;„P®"], 



with Xi,yi > 0. The relevant hyperplane is given by 

iv[p.«(«+(i-^)p)%]=E(:)(i-f)'..=o. 

1=1 ^ ^ 

Renormalising Xi, this can be rewritten as (")(— l)*ij = with ij = (| — lYxi. Next we will show that this 

hyperplane touches the set of separable states by constructing a set of 2n+ 1 separable states spanning the hyperplane. 
The first n+1 states are obtained by twirling 

iV')fc = |ol)|oo)®'=|ol)®"-^ 
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for k = 0, - ■ ■ ,n. The twirled state will be !7L'^*(T^T^*F)-invariant and will satisfy Xi = for all i and yj ^ for j < k 
and ijj = for j > k. The last n states are obtained by twirling 

1 

1^-)^= E l^J)NJ)loo)«"-^-l|ol)«^ 

for fc = 0, • • • , n — 1. The twirled state will be C/J7*(V^iV^*F)-invariant and the x coordinates will satisfy Xj = for 
j > n — k and Xj ^ for j < n — k. Therefore the coordinates of the 2n + I states are linearly independent and 
the choice of \ip)k guarantees that the states will lie in the hyperplane. An interior point in the convex hull of these 
2n + 1 points can be obtained by choosing xq = e and yj = Xj ~ e, for e sufficiently small. One verifies that this 
point belongs to the hyperplane using the identity Y^"^q ~ From this follows that the conjecture is 

equivalent to the statement that the J7[/*(t/iy;*i^)-invariant states with coordinates yi — e , Xi — e and xi — e + 5, 
where j = 1, • • • , n and i = 0, 2, • • • , n, are entangled for e > small enough and all 5 > 0. 



C. Discussion 

The general separability problem has been proven to be NP-hard [2^. However, above we showed that the dis- 
tillability problem can be reformulated as the question of entanglement of a particular set of one parameter states. 
The first and by far still the most elegant tool for detecting entanglement is the partial transposition criterion |29l |. 
However, it is not useful to solve the dual entanglement problem, as it was proven in j^l that W^^n is a decomposable 
entanglement witness if and only if p has positive partial transition. Hence W tgn will be either a non-decomposable 
witness or no entanglement witness at all depending on whether is 1-undistillable or 1-distillable. Thus in order 
to solve the dual problem, we will either need some powerful tool for detecting PPT- entanglement or a tool for proving 
separability. However, in the latter case, the original formulation stemming from Theorem II. II seems easier and for 
this purpose we present an efficient algorithm for detecting distillability in section Hvl 

Probably the most powerful method for detecting entanglement is the complete family of separability criteria 
introduced by Doherty et. al. [s^jl^- Basically, their method relies on a hierarchical characterisation of separable 
states which they use to devise a computational algorithm for detecting entanglement. It has a number of very 
appealing features: (i) The set of criteria is complete, all entangled states are detected at some stage, (ii) The 
criteria can be cast into a semidefinite program which is a convex optimisation problem for which efficient algorithms 
exist, (iii) When a state is found entangled, their algorithm automatically yields an entanglement witness for that 
state. These entanglement witnesses turn out to be of a special form, namely such that after some manipulation 
the associated bihermitian form can be written as sums of squares. Bihermitian forms which can be written as sums 
of squares are canonical entanglement witnesses. Therefore it is always possible to extract some analytical provable 
entanglement witness from the output of the algorithm. In principle therefore, the distillability problem can be solved 
for any number of copies using our dual formulation together with the algorithm associated with the complete family 
of separability criteria. 

Another numerical method based on semidefinite programming is the one introduced by Eisert et. al. |32l |. There 
the separability problem is cast as a global optimisation problem with polynomial constraints. Using the theory of 
semi-definite relaxations, a hierarchy of efficiently solvable approximations to the optimal solution is provided. Every 
entangled state is necessarily detected in some step of the hierarchy and since a global optimum will be found at 
some point, separability can also be detected. However, it should be noted that the reformulation of the distillability 
problem as a separability problem is unnecessary here, as their method allows for a simple way of checking whether 
or not an operator is an entanglement witness. Using analogous techniques, it can be extended to testing whether or 
not an operator is a Schmidt number 3 witness [b^ . 

The family of criteria introduced by Doherty et. al., however, has the advantage that an analytical canonical 
entanglement witness can be extracted. Unfortunately, we were unable to test either algorithm, as both methods 
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FIG. 4: The probability of detecting 1-distillability as a function of the test number for random states drawn from T) xlA. For 
visual purposes only the first 100 tests are depicted and the points have been smoothened out to curves. 



seem only practical for low dimensional systems. The interesting case, two copies of the Werner states for d — i could 
hence not be tested. It would be worth investigating whether these numerical methods could be simplified for the 
highly symmetrical states we are interested in. 

Numerical solutions are one option, another possible approach would be one of a more indirect nature. A powerful 
method for proving that a certain state is entangled is to show that, when shared by two parties, the state can 
enhance typical quantum operations such as tclcportation or distillation. In particular, it was shown that PPT- 
preserving operations can be simulated using LOCC operations when both parties share PPT entangled states [2o| . 
The class of PPT-preservin g op erations is strictly larger than the LOCC class, so that we can expect it to do more. 
Recently it has been shown |2l| that every entangled state can enhance the so called conclusive teleportation fidelity 
of some other state. It was also proven that for every PPT state a there is a 1-undistillable state p such that p® a is 
1-distillable. These characterisations of entangled states seem very promising as a way of proving that a certain state 
is PPT entangled. In Ref. [3 ^ class of PPT states was constructed which was shown to provide overall convertibility 
of pure entangled states. In particular it was shown to be able to increase the Schmidt number of a pure state. Now 
the PPT entangled states we obtained for one copy of the Werner states (tetrahedron ABCO in Fig.|3Jl are of a similar 
form, and although not necessarily in an optimal way, they too can help increase the Schmidt number of a pure state. 
Similar activation effects can be expected from the conjectured PPT states derived from two copies of the Werner 
states. 



IV. THE PEASANT'S METHOD 

In this section we will outline a powerful algorithm for the detection of distillable states. As n-distillability of p is 
equivalent to 1-distillability of p®", it is clear that we can confine ourselves to the study of 1-distillability. The basic 
problem we need to solve is the following minimisation over Schmidt rank two vectors: 

min ('(/'Ip^^I?/') < 0. (5) 

The numerical method employed by Diir et al Il3 | converts this problem to a minimisation of the minimum eigenvalue 
of certain matrices. Their method involves also calculating an inverse of a matrix in every step. The peasant's method 
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[48| does not require this and is also conceptually simpler. 
We can rewrite equation ^ as 

min Ti ( \ih) ( ihl p'^'' ) = min Tr(Dp^«). 

It is clear that the operators D = pl^'Y the role of a distillability witness, in much the same way entanglement 

witnesses detect entanglement j^^. So unlike in the entanglement problem, we have a complete characterisation of all 
distillability witnesses. Now it is well known that the positive map associated with an entanglement witness detects 
much more entanglement than the witness itself js^. The positive maps D(-) associated with D = |'0)(V'l can be 
chosen to act as 

D{p) = 

with P ~ \0){a\ + |1)(6|, with (6|a) = 0. This in effect is a generalisation of Theorem lI.il stemming from the fact that 
in n (8) 2 all NPT states are distillable (see also 0,0,E3)- Thus the problem is reduced to checking whether there 
exist vectors \a), \b) such that a has a negative eigenvalue. One way of doing this is to parametrize a countable subset 
of vectors which is dense within all vectors, such as the one introduced in Ref. Explicitly their set takes the form 
G = {E A.N)|(Ai, A2, • • • , Xd) e Gn} with 



Gat = < — e =1,,— e , ■ ■ ■ , \ - 1^ ^-i ] } 

ior < Pi < Qi < N;0 < Vi < Si < N. Thus for every N, we can construct sets of pairs \a) and \b), and taking N 
increasingly large we will detect all 1-distillable states, except those arbitrary close to the boundary of the convex set 
of 1-undistillable states. 

In practice two improvements can be made which greatly enhance the performance. First note that the countable 
subset above will yield vectors \a) and \b) not necessarily orthogonal. Furthermore, it is clear that such a countable 
subset will in general not pick vectors uniformly distributed according to the Haar measure. One way of overcoming 
this is to take for \a) and |6) two columns of a random unitary. Here orthogonality and uniformity are automatically 
guaranteed. The associated algorithm works very well for low dimensional density matrices. For higher dimensions, 
a local optimisation of the \a) and \b) yielding the minimum after a certain cut-off number of tests turns out to work 
well. 

We have checked the distillability of the f/t/V^Fi^-invariant states over the complete range of parameters for 1 and 
2 copies for d = 3. We easily recovered the proposed boundaries for 1-distillability. For two copies, the states act on a 
C(;4 (g)C(ii Hilbert space and numerical matrix manipulations in a space of this magnitude seem very hard. Fortunately, 
the states are very sparse and the peasant's method only requires minimisation of the minimum eigenvalue of a 2(P x 2(P 
matrix. We were able to detect distillability for the Watrous states in e > 1/d — 1/2 readily and exhaustive testing 
suggest that also the proposed boundary for 2-distillability is correct. This provides very strong evidence that the 
Werner states are 4-undistillable for /? > —1/2. 



A. Volume of 1-distillable states 

As an application of the peasant's method we will give a numerical estimate of the volume of 1-distillable states 
for low dimensional quantum states. A similar numerical estimate has been carried out for entangled states 0,0. 
When talking about volumes on the set of density operators, it is clear that the results will depend on the measure. 
We choose the measure applied in Ref. [i^ as it seems very natural. 

A general bipartite quantum state p acting on Ha ^'Hb can be expanded by virtue of the spectral decomposition 

as 



p=UDU\ 
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FIG. 5: Probability of a random NPT state drawn from V x U to he 1-undistillable. The curve has been drawn to guide the 
eye. 

with D = (da) diagonal and U unitary. The measure we are going to use is the product measure x U. Here V 
represents the uniform distribution of the points on the manifold given by da = 1. A simple method for generating 
such a distribution from independent uniformly distributed random numbers chosen in the interval (0, 1) is outlined 
in Appendix A of Ref. (3ij |. Similarly, U is chosen to be the uniform measure on imitary matrices (the Haar measure). 
To generate random unitaries according to this measure one can use the algorithm from Ref. |4l[ |4^ which relies on 
a decomposition of a general unitary in two-dimensional unitary transformations. A simpler method for generating 
Haar unitary matrices is as follows Take a random matrix A, whose entries are complex numbers which 

are independently normally distributed with mean zero. The QR factorisation A = QR, such that R has positive 
elements on the diagonal, then yields Q distributed according to U. 

We have tested the peasant's method on 10^ density matrices acting on C' ® for d = 3,4,5,6,7. A striking 
feature is that the vast majority of the distillable states were detected in the first few tests. For example for d = 3 
about half of the NPT states are found to be distillable in the first test. The estimated probability of success as a 
function of the test number is displayed in Fig. 0] For d = 7 the probability of finding a distillable state in the first 
test is about 1/6. This seems to suggest that the volume of 1-distillable states drops to zero for high dimensions. 
The opposite turns out to be more likely. In Fig. O the probability of an NPT state being 1-undistillable is plotted 
versus the dimension d. To obtain sufficient precision we carried out 10^ random tests per state, and in addition 
lO'^d optimisation steps seeking for a local minimum. Of course, this method does not guarantee to detect every 
1-distillable state, but we obtain an upper bound of the number of undistillable states. Note the distinct peak at 
d = 3, 4, the reason for this will be explained elsewhere 0. In Fig.|Hlthe same graph is drawn, but now PPT states 
are included. 

It is tempting to conclude from this numerical evidence that bound entangled states are primarily a phenomenon 
present in low dimensional quantum systems. In high dimensional systems most undistillable states are therefore 
situated in the immediate neighbourhood of the set of separable states. Our results are consistent with the fact that 
bound entanglement for continuous variables is a rare phenomenon [4^. In particular it was shown that the subset 
of undistillable states is nowhere dense in the set of bipartite continuous variable states. From this it follows that the 
set of undistillable continuous variables states does not contain any open ball, an argument which was made explicit 
for separable states in 0. Given an infinite dimensional separable state one can construct sequences of closer and 
closer states all of which are entangled. Following the same methods one can explicitly construct a distillable state 
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FIG. 6: Probability of a random state drawn from V x U to he 1-undistillable. 

in any e-neighbourhood in the trace norm of any state 33]. However, note that in Ref. a parametrized family of 
measures on states was introduced, which in some region yields states primarily 1-undistillable. 

V. CONCLUSION 

The main result of this paper was that the distillability problem can be formulated as just a special instance of the 
separability problem. We have discussed several ways in tackling this separability problem which we believe merit 
further study. We have outlined an efficient numerical method for detecting distillability, and provided strong evidence 
that the distillability conjecture is valid at least up to 4 copies of the Werner states. The method was also used to 
make an estimate of the volume of 1-distillable states for = 3, . . . , 7. 

APPENDIX A: THE RAINBOW STATES 

Let us consider (we omit the indices {1, 2} and {3,4}) 

de — 1 „ me — 1 „ 1 — (m + d)e + dmS „ „ 

p = 1„ ® Id +—J— Irn (»Fd + Fra «) Id + ^ (g) Fd- 

d m dm 

In what follows we will assume that 3 < m < d. The set of density operators is restricted by the following inequalities 

md md 
1 ^ m+d ^ 1 l^p 
md 7nd m d ~ 

md md 

The partial transpose is given by 

= Qm <E)Qd + meP„i ®Qd + deQm ® Pd + mdSP„i (g) Pd, 
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and it is easy to see that the states are NPT iff e < or 5 < 0. Now we wih show that these states also include some 
PPT entangled ones. From Lemma [ill. 1 1 and the results in liS] it follows that 

mP„ (g) (Id- — Pd) 
m 

is an entanglement witness for p-^^ . The partial transpose of an entanglement witness is again an entanglement 
witness: 

F^(g){ld- — Fd). 
m 

So that in addition apart from the states e < or (5 < 0, the states that satisfy 

err?{d^ - 1) + dm5{m - d) < 

are entangled. 

Let us now look at the distillation properties. Similar arguments to the ones we used for the UUVV F-uiYa,Tia,nt 
states apply here. For e or J sufficiently small, one can find n-undistillable states. The following vectors provide 
apparently the optimal boundaries for the 1-distillable states (we will only be interested in the states having S > 0). 

1. \^) = |00)A|00)i3 + |10)a|10)s gives rise to 

2 + 2{de - l)/d + A{me~ 1)/to + 4(1 - {m + d)e + dm5) / {md) < 

2. \^) = |00)A|01)i3 + |10)a|11)s gives rise to 

^<^-i (Al) 



m 2 



Now let us take two pairs 



Pi - Jl-m ® '^d H 1 H ^ni ® "I 1 ® : 

a m dm 



and 



78^1,56 , de-l ^78^c.56 , ^^^^ " ^ p78 ^ ^ 56 , I - {m + d)e + dm5 ^ ^gg 

m 

1,5 ^ n2,6 



Taking both pairs together, and projecting upon P^' (g) P^' , we end up with a C/C/FV^P-invariant state in m m 
with 



^, _ e{d^5 + d^-l) 
5' ^ 



^2g2 + rf2 _ 1 
e2(d2_l) + rf2^2 



^2g2 + rf2 _ 1 • 

Now we know when these states can be distilled. The states distillable with this protocol are shown in Fig. [71 One 
verifies that projecting upon P^f ® P^^f performs worse. The set 1 — ^ + e-^^^ ^^^m^2^^ contains states 
of all kinds: 1-distillable, n-undistillable but n + 1-distillable, NPT undistillable (conjectured), PPT bound entangled 
and separable. 
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